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ABSTRACT 



Antarctica temperature variations have long puzzled 
researchers, giving them reason to believe that a mean 
diurnal temperature variation existed during the totally 
dark polar night. It was later found that diurij^l varia- 
tions were most evident when the only days averaged were 
those having low wind speeds and little cloud cover. This 
led to a conclusion that a physical explanation was not 
necessary and that the variation could be explained as re- 
sulting from a statistical bias. 

This study is a continuation of an investigation begun 
at the U. S. Naval Postgraduate School in 1963 to determine 
whether or not a true Antarctica diurnal temperature varia- 
tion exists and to determine the extent of the statistical 
bias effect in modifying the mean temperature variations of 
specially selected days. In broadening the original study, 
the authors utilized an additional 12 years of data, a se- 
cond geographical location and more accurate methods of 
obtaining results. 
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1. Introduction. 



Many investigators, including Simpson, Court, Rouch and 
Hisdal, have observed an apparent diurnal temperature varia- 
tion during the totally dark polar night season at various 
stations in Antarctica. In general, the averaged daily tem- 
perature curves exhibited minimum temperatures near noon, 
local time, and maximum temperatures near midnight. Intui- 
tively, a diurnal temperature variation is unexpected in a 
region which is receiving no insolation and much effort was 
expended in finding a physical explanation for this phenome- 
non. 

Simpson spent considerable time investigating the diur- 
nal temperature variation he found at McMurdo Sound in one 
year’s data from a Scott expedition. [_ll] However, he wrote 
in 1919 that he could report only negative results in his 
attempt to explain the strange temperature behavior. 

Rouch, in 1941, could not find a satisfactory explanation 
for the daily temperature variation he found in polar night 
data from both Charcot and Byrd expeditions. 

Other investigators have, more specifically associated 
the diurnal temperature variation with certain types of days. 
Namely, they found that days with clear skies and light or 
moderate winds had mimimum temperatures near noon, local time 
whereas, overcast and windy days had a maximum temperature 
near noon. However, Court in 1951, was still unable to ex- 
plain the phenomenon. Q> 3 

In a study of the effect of the moon on atmospheric pres- 
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sure variations, Bartels found a similar and also unexpected 
diurnal pressure variation associated with specially selected 
days. [2]] As related in his unpublished Gottingen disserta- 
tion in 1922, he determined that the pressure variation phe- 
nomenon was the result of a statistical bias induced by the 
day selection process. Bartels’ statistical bias theory re- 
ceived additional support from van der Bijl as a result of 
his investigation of the diurnal pressure variation in north- 
western Europe during specially selected days. &2] 

It remained, however, for Hisdal in 1960, to independ- 
ently ’’rediscover” that the statistical bias effect could also 
explain the polar night diurnal temperature variations associ- 
ated with the above mentioned specific types of days. D1 

Because of the extremely cold surface in Antarctica, 
periods with light or moderate winds, and hence little turbu- 
lence are generally characterized by a temperature inversion. 
Similarly, periods with little or no cloudiness, and hence 
with little or no downward radiative warming, are also gener- 
ally characterized by a temperature inversion. For the oppo- 
site reasons, non- inversion periods can be associated with 
windy and/or cloudy days. Periods during which a temperature 
inversion exists will thus usually exhibit lower temperature 
averages than non- inversion periods, but more importantly, 
temperature records for such periods will have similar curva- 
ture . 

Figure 1 shows a hypothetical temperature curve for a 
period of seven days. Included are two periods characterized 
by temperature inversion. These periods are of unequal length 
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Figure l; Hypothetical temperature curve for seven 
days showing periods of inversion and non- inversion. 
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TIME 





and they do not coincide exactly with any 0000 to 24000 day. 
They do, however, have lower average temperatures than the 
non-inversion periods and their lowest temperatures occur 
near the middle and their highest temperatures occur at or 
near the end of the periods. In general, the temperature 
curves during the inversion periods are concave. 

If all 24-hour periods, having given starting times and 
each totally within a period of inversion are selected from 
the record and the resulting mean daily (24-hour) temperature 
curve computed, it will be found to have the general concave 
shape of the longer inversion periods. If the starting time 
of the selected 24-hour periods is shifted, the time of occur- 
rence of the minimum temperature of the computed mean curve 
will also be shifted. The shape of the mean daily temperature 
curve is thus partially determined by a statistical bias in- 
duced by selecting only specific 24-hour periods which indi- 
vidually tend to have like curvature. 

For inversion-selected days, the mean temperature varia- 
tion curve is therefore, composed of both the meteorological 
temperature variation, if any, and a statistical bias. His- 
dal felt that this bias could be of sufficient magnitude to 
mask the true meteorological temperature variation. By im- 
plication, a sufficient number of inversion type days over a 
long period could result in an apparent mean diurnal varia- 
tion in investigations where all days are considered. 

Barrigar commenced a study at the United States Naval 
Postgraduate School in 1963 to determine the magnitudes of 
the meteorological temperature and statistical bias variation 
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of the mean temperature curves computed for specially selected 
days using six years of polar night data from McMurdo Sound, 
Antarctica. In addition* he computed the mean temperature 

variation using the data for all days during the totally dark 



This study is a continuation and expansion of that initi- 
ated by Barrigar. It utilizes another four years* data from 
McMurdo Sound and investigates a second geographical location 
with eight years* data from the South Pole. In this study, 
an additional method is employed to obtain a more accurate 
representation of the statistical bias than was obtained pre- 
viously. A more thorough investigation is also made of the 
applicability of the statistical bias effect over time in- 
tervals greater and less than 24 hours. The authors briefly 
investigated the varying durations of the inversion periods. 
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2. Data and Processing. 



The data, in the form of punched cards, used in this 
investigation was obtained from the National Weather Records 
Center, Asheville, North Carolina. Each punched card contains 
the meteorological information of one surface observation in 
the format described in the "Card Deck Reference Manual, 144 
WBAN 1". Cards for the periods January, 1962 through October, 
1965 for McMurdo Sound (78° South, 166° East) and January, 

1957 through January, 1965, for Amundsen Scott station at the 
South Pole were received. (Figure 2 shows the locations of 
these stations.) Most of the card-recorded observations 
were at three-hour intervals although during some periods, 
hourly observations were available. Some observations were 
completely missing from Amundsen Scott station. This was 
particularly true in a random fashion during 1957 and again 
in 1964 when many of the 1500 observations were missing. In 
addition to this data, a magnetic tape from the study con- 
ducted by Barrigar was available. This tape contained data 
from McMurdo Sound at three-hour intervals during the polar 
night period, 15 May through 15 August, for the years 1956 
through 1961. 

Solar altitude curves prepared by Laevastu indicated 
that the duration of the dark season at McMurdo Sound is 
from 20 March through 20 September. (sQ However, to eliminate 
twilight effects and any incoming solar radiation due to un- 
usual atmospheic refraction conditions, the shorter periodT in- 
dicated above was used to insure total darkness at stations 
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during the period of investigation. 

A Fortran-63 program was written for use with the Con- 
trol Data Corporation 1604 (CDC 1604) digital computer to 
transfer the data from the McMurdo Sound tape and the new 
data cards .to a new magnetic tape containing ten years of 
polar night data at that station. 

Although this investigation of the Amundsen Scott 
(South Pole) data examines just the same three-month polar 
night period as that investigated at McMurdo Sound, all of 
the data at three-hour intervals for the South Pole was put 
on magnetic tape. This will enable future investigators to 




Geographical location of McMurdo Sound, Antarctica. 
Amundsen Scott station is located at the South Pole. 
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examine the South Pole data for any time interval. 

The Fortran-63 program which was written to transfer 
the South Pole cards to magnetic tape incorporated several 
selective features. To simplify further processing, it 
transferred only observations at three-hour intervals, and 
it did not transfer cards which had missing temperature, 
wind, or total cloud amount data. This eliminated about 
100 observations. In addition, the program converted all 
temperatures from three to four columns. This was necessary 
to make temperatures less than minus one hundred degrees 
compatible with the other recorded temperatures. 

The magnetic data tapes were subsquently processed on 
the CDC 1604 digital computer, using Symbolic Coded Relocat- 
able Assembly Program (SCRAP) to: 

(1) compute for each 24-hour interval over 
any desired period (e.g. 15 May to 15 
August) the mean values of temperature, 
cloudiness, atmospheric pressure, and 
windspeed for each observation time 
within the 24-hour interval; 

(2) compute the total mean values determined 
in (1) for all the years processed; 

(3) compute the means in (1) and (2) for 
intervals of less than or greater than 
24 hours; 

(4) compute the means in (1), (2), and (3) 

for various starting times of the interval; 

(5) select the intervals for which the above 
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computations are made on the basis 
of whether or not all of the ob- 
servations in the interval met the 
specified windspeed and cloud amount 
criteria ; 

(6) and compute the Lamont corrected values, 
described in section 3, of the means 
above . 

The data tapes were also processed using Fortran-63 and 
SCRAP programs to compute a mean daily temperature for each 
specific month and day averaged over all the years possessed. 

All other computations, such as column and diagonal 
averaging procedures detailed in section 5 were done by cal- 
culator . 
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3. The Lamont Correction. 



The Lamont correction is a technique for modifying the 
values of a continuous parameter so that the values are equal 
at the beginning and end of a specified time period. This 
technique simplifies the study of cyclic variations by remov- 
ing the effect of long term trends from the short term cyclic 
variations. The Lamont correction was utilized in this tem- 
perature study to remove the seasonal trend from the diurnal 
cycle . 

Normally the correction is applied by modifying all val- 
ues of the parameter about the value at the center of the cy- 
clic time period to equalize the end values. This produces 
the graphic effect of rotating a curve about its central val- 
ue until the end values are equal. The process alters the 
absolute values of the parameter without affecting the arith- 
metic mean of the time period provided the long term trend is 
linear. A more complete discussion of the Lamont correction 
is contained in Conrad and Poliak. £ 4 ] 

The Lamont correction applied in the computer sub-routine 
mentioned in section 2 modifies all values during the cyclic 
period by an amount proportionate to the difference between 
beginning and end values of the parameter. This effectively 
rotates the curve about the starting value until the end 
value equals the starting value. This deviation from the nor- 
mal procedure results in the alteration of the absolute val- 
ues of the parameter and a change in the arithmetic mean. 
However, relative values of the parameter are preserved as 
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before and since only relative values of the cyclic variation 
are required, this simpler method is utilized. 

The Lamont corrected temperatures were computed using 
the following equations: 



T iL 

%<■ 




(l) 

co 



where: i = sequential numbering of the observation 

times in the cyclic period (i=0,l,2, . . . . ,n) 

n = number of observation times in the cyclic 
period 

= observed temperature at i th observation 
time 

Tib = Lamont corrected temperature at i th 
observation time 

g^ = Lamont correction to be applied at i** 1 
observation time 

As an example of Lamont's correction procedure, appli- 
cation of equations 1 and 2 to the hypothetical observed tem- 
peratures in table 1 yields the g^ and T^l values tabulated. 



TABLE 1 



Hypothetical temperature observations (T^) 
ing Lamont corrected temperatures CT^)- 


and 


the result 


- 


Hour 


of Observation 


00 


06 


12 


18 


24 


i 
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1 


2 


3 


4 


T i 




2 


3 
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12 


Si 
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-2 


-4 
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T iL 




2 


1 


2 


3 


2 
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A graphical representation of the Laraont process utili- 
zing the information in table 1 is shown in Figure 3. The 
observed temperature curve with a long term upward trend and 
Lamont corrected temperature curves which pivoted about both 
the central and initial temperature values are plotted. Note 
that the Lamont curves have different absolute temperatures 
but identical shapes. 
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Figure 3 



Hypothetical observed and Lamont 
corrected temperature curves. 
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4. Theoretical Discussion of the Separation Technique. 

The purpose of this section is to examine a method for 
separating the mean temperature variation record into its 
meteorological and statistical components. 

The apparent diurnal temperature variation obtained for 
special types of days during the polar night may be considered 
to be the sum of the daily meteorological variation and a sta- 
tistical bias. Since a purpose of this study is to determine 
an estimate of the variability of these two components, it is 
necessary to examine how the separation may be performed. 

Consider first the composition of the observed tempera- 
tures. Table 2 tabulates a hypothetical statistical bias for 
each six-hourly observation time in a 24-hour period. 

TABLE 2 

Hypothetical statistical bias for each six-hourly obser- 
vation time in a 24-hour period. 



OBSERVATION TIME i 
0 th 



TEMPERATURE 

3 



1 



st 



1 



2 nd 



0 



3 rd 



1 



4 th 



3 



As pointed out by Hisdal, the statistical bias is not a 
fixed value for each like observation time. Instead, it de- 
pends on the starting time of the special type of 24-hour 
period selected. If the selected day, for instance an inver- 



28 



sion day, begins at 1200, the statistical bias will be in 
phase with the periods' starting time and in the example of 
table 2, the 1200 statistical bias will be three degrees. 

The statistical bias is thus comparable to a free wave which 
moves through the meteorological temperature curve with its 
influence on the observed temperature record, dependent upon 
the starting time of the period selected. 

The actual composition of the total observed temperature 
curves for days (24-hour periods) with various starting times 
selected on the basis of their being characterized by an in- 
version is illustrated in table 3, and in Figure 4. In table 
3, the hypothetical Lamont corrected meteorological tempera- 
tures (Tjp) of table 1 are added to the statistical bias 
temperatures of table 2 at the four starting times to pro- 
duce the total observed temperatures. Utilizing the values 
in table 3, Figure 4 illustrates the procedure graphically. 
The upper curve is the meteorological temperature curve 
which has a constant daily variation. Below it are plotted 
the statistical bias curves and the resulting observed tem- 
perature curves at each of the four starting times. It is 
readily seen that the observed temperature curves are no 
longer representative of the meteorological temperatures and 
that the time of minimum observed temperature shifts as the 
statistical influences are shifted due to the different 
starting times of the special days that were selected. 

Now that the composition of the observed temperature 
records has been examined for specially selected days over 
various 24-hour periods, the next step is to examine how the 
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TABLE 3 



Observed temperature variations resulting from the meteoro- 
logical and statistical bias variations from tables 1 and 2 
respectively, for various starting times. 



Starting 


time: 


0000 






Hour : 
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Meteorological : 


2 


1 


2 


3 


2 


Observed : 


5 


2 


2 


4 


5 


Starti 


ng 


time : 


0600 






Hour ; 


06 


12 


18 


00 


06 


Statistical : 


3 


1 


0 


1 


3 


Meteorological : 


1 


2 


3 


2 


1 


Observed : 


4 


3 


3 


3 


4 


Starting 


time : 


1200 






Hour : 


12 


18 


00 


06 


12 


Statistical : 


3 


1 


0 


1 


3 


Meteorological 


2 


3 


2 


1 


2 


Observed : 


5 


4 


2 


2 


5 


Starting 


t ime : 


1800 






Hour ; 


18 


00 


06 


12 


18 


Statistical : 


3 


1 


0 


1 


3 


Meteorological: 


3 


2 


1 


2 


3 


Observed : 


6 


3 


1 


3 


6 
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I 

Figure 4: Total observed temperature variation curves 
for various starting times constructed from the given 
meteorological and statistical bias variation curves. 
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original values of meteorological variation and statistical 
influence may be extracted from the observed temperature re- 
cord. To facilitate the explanation, the previously used 
hypothetical data has been put in the form of table 4. The 
values are recorded in the following form: 

t^; where t = hypothetical observed temperature 
s = statistical influence 
m = meteorological influence 



TABLE 4 



Hypothetical observed temperatures 
starting at 0000, 0600, 1200, and 


for various 24- 
1800. 


-hour 


periods 


00 


06 


12 18 

Observation 


24/00 

Time 


06 


12 


18 


T^x <3 ^ 
E 




^ 22 " 




<2 








M , 
















P 3 1 

E ^2 

R 






\3° 3 \3j - 








A -0 

T ^2 

TJ 


2*1 




V 4 1 


2° ^ 
\J2 




5 2 

d 3 




R 3 1 

E ^2 


1? 


3 2 X 


\ .3 \ 

\ 3 


ol \ 
v 3 2 x 




O 1 

\ S 2 


x a 3 
X3 








X, 


\ 


\ 


^ 


\ 



The separation technique is that used by van der Bijl 
in his study of diurnal pressure variations. It consists 
of averaging the values of table 4 in two ways: by columns 

to determine the meteorological influence and diagonally to 
determine the statistical influence. Once either of these 
two averages has been computed, the other may be determined 
by subtraction from the average of the corresponding observed 
values . 

As detailed by Barrigar, these methods are shown below. 
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Notation is as follows: 



trk,x where t 



temperature (F) 



r = (when present) denotes relative 
temperature 



k,x = m,hh refers to meteorological 

influence at the indicated hour 
of observation 



k,x = s,i refers to statistical in- 
fluence at i^h observation 
time in the period 

hh = hour of observation. 



Averaging the observed temperatures by columns, and sub- 
tracting the constant values of the statistical bias, results 
in the meteorological variation. 





Averaging the observed temperatures and subtracting the 
constant value of the meteorological influence, results in 
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the statistical influence. 



s.o ~ 






_ 2.0 

H 


_ 8 
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r 3 , 


i,l ' 






_ /* 

H 
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V 


= / 


s.i s 


{l±l±i±L} 


- e^> 


„ 8 

H 


8 


= 0 


VJ = 




- 


_ 12 


8 


= / 








- lo 

H 


8 


* 3 



In this example, the originally assumed values of the 
meteorological and statistical influences have been recovered. 
In actual practice, however, the only known values are those 
of the observed temperatures. The constant terms just sub- 
tracted from the observed temperature values are unknown but 
they may be determined as follows: 

Using the additional notation: 

C m = a constant which is dependent upon 
the meteorological influence 

C g = a constant which is dependent upon 
statistical bias 

the previously computed averages become: 

'3 



-6 



m,oo 
' 1 * 1,00 



~ C 



j- - r 



H 






-%■ - c, 






i 






H 

12 

H 

*/ 



±L - r 

U V- < 



C. 



(-S 
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Using 0000 as the starting time, the previously com- 



puted diagonal averages become: 
20 



■*s,o 

K* 



TT - C. 



c - 

YV% 



H 

- ^ - c 

- 4 - c 

- 11. - f 

H U 

Zo 
4 



Vv) 



Zo - C 






W| 



►vi 



W\ 



yv> ,oo 






muz 



+ 4 



*V8 



Making the appropriate substitution for C m in the last 

set of equations: 

/ — - 4 
^5,0 " 4 Lj 

= -J- 17-il +iL+.Q+iIN - url 

H V [A H H H H J 7 L SJ 

r 12 ^ ( ^S.n + ^,,+ j&j 

c IT ~ (^ 5 -° + ^s,i + 

Similarity, for the starting times at 0600, 1200 and 
1800, we obtain, respectively: 



^S,0 - 3 K\ * ^ **^$,3 3 1 



^s,o + ^s,/ ' 3 ^ 5j z +^ S)3 = 

^S,0 r ^s,l + ^S,z " 



= 5" 



= / 



0) 

& 
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Equations ( 3 ) through (6) comprise a set of four linear, 
homogeneous equations in four unknowns, which can be solved 
using the usual mathematical techniques. The solution yields 
the following statistical bias temperatures: 

t s , 0 = 3 

V 1 = 1 

t g , 2 = 0 

t s , 3 - 1 

t s , 4 - ts , 0 = 3 . 

These are the originally assumed values and subtracting 
from the appropriate observed temperatures yields the origi- 
nal meteorological temperature. 

By solving sets of equations, it is thus possible to 
determine the actual values of the statistical and meteoro- 
logical constituents. However, it will be shown in the next 
section that for the purposes of this study, it is not neces- 
sary to determine the constants C s and C m . 
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5. Discussion of the Separation Technique Used in This Study 



The purpose of this section is to illustrate the techni- 
que employed in this study for the separation of the mean ob- 
served temperature variation curves for selected periods into 
their meteorological and statistical components. The SCRAP 
program mentioned in section 2 computes the average observed 
temperature at each observation time for the "days" selected 
of specified length and starting times. These are represented 
by the observed temperatures listed in table 4. The hourly 
meteorological and statistical contributions are unknown. 

As shown in the previous section, it is necessary to solve a 
set of linear homogeneous equations to determine the specific 
values of statistical bias and meteorological effect. In 
fact, the equation sets would consist of eight equations in- 
stead of four as illustrated in the example. This is a re- 
sult of our determining the mean temperatures at three-hour 
intervals instead of the examples' six-hour intervals in order 
to obtain more detail temperature variation curves. 

As this study is concerned with the shapes of the tem- 
perature curves and not with the absolute temperature curves 
and not with the absolute temperature values, it is suffi- 
cient to determine the relative values of the meteorological 
and statistical temperature constituents. 

By using the procedure of section 4, the observed tem- 
peratures in table 4 are column averaged to determine the 
meteorological influence. If the constant C s values are now 
neglected, the remaining numerical values become relative 
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values of the meteorological temperature at each hour as 
shown below where the subscript r denotes relative values. 



i 

-b 

-L 



Column Averages Relative Meteorological 

Values 



rr\ ( 00 


13 

- H 


-C t 


^rmoo 


13 
- 4 


- 3.2 5 




_9_ 
" H 


- C s 




<r| J- 

II 


- 2.25 




_I3_ 
r H 


- C s 




/3 

- H 


= 3.25 


'Vfi 


17 

- H 


-C s 




ll 


= 4. 75 




_/3 

= H 


-C s 


4 


H 

- H 


= 3.25 



Subtracting these relative values from the appropriate 
observed temperatures, yields the relative value of statisti- 
cal bias for each hour as shown: 

+ n>0 --S -2.25- L2S 

^ v = 2-2.2r, -.25 

i n/l = 2 - 3.2 S = -I.2S 

-L n l -- H - 4.25 = -.25 

ij - 5 ~ 2, 2 iT - 7*5” 
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The originally assumed and just determine meteorologi- 
cal and statistical bias temperatures are plotted in figure 
5. Although the numerical values of the original and rela- 
tive values are different, the shapes of the curves are 
identical. The above example utilized the observed tempera- 
tures of the period starting at 0000, but the method is appli- 
cable for any period. 

Note that in this procedure the relative meteorological 
temperature was determined before computing the relative 
statistical bias. It can be shown that by averaging diagonal- 
ly, eliminating the constant C m and making appropriate sub- 
tractions, the relative statistical bias can be determined be- 
fore the meteorological temperature. 

Because of the existence of random errors in the data 
utilized in the study, it is desirable to use computations 
which minimize their effect. Barrigar showed that the best 
method of obtaining the meteorological influence is by col- 
umn averaging and that the best method for obtaining the sta- 
tistical bias is by the diagonal averaging method. In fact, 
using the diagonal instead of the column averaging procedure 
reduces the variance in the computed statistical bias by a 
factor of seven where the observations are three hours apart. 
Both averaging methods are used in this study. 



39 



Figure 5 



Assumed and computed relative values of 
meteorological and statistical bias tem- 
peratures. 
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6. General Discussion of the Results. 



The results for the McMurdo Sound and South Pole sta- 
tions are individually presented in the following two sec- 
tions. As described previously, the data reduction and pro- 
cessing led to a series of temperature variation curves for 
each of the stations. A set of curves consisting of the 
mean observed meteorological and statistical bias tempera- 
tures for the specified intervals has been drawn for each of 
the various selection criteria utilized. Each of the sets 
of curves consists of a Figure a and a Figure b. The meteor- 
ological and statistical bias variations of the a Figures 
were computed using the previously mentioned column averag- 
ing method. Figure b variation curves were computed using 
the diagonal averaging method. The a Figures thus contain 
the most representative meteorological temperature curves 
and the b Figures the most representative statistical bias 
curves. In each instance, the mean temperature variation 
curve is that of the interval beginning at 0000. 

The frequent significant variations in the curves com- 
puted using the two methods justifies the additional diagonal 
averaging computations which were not utilized in the pre- 
vious study. As an example, note the difference in the sta- 
tistical bias curves of Figures 17a and 17b. 

For easy reference, tables 5 and 6 list the number of 
periods meeting the various selection criteria at the MuMurdo 
Sound and South Pole stations, respectively. In comparing 
the curves which follow, it is important to realize that the 
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same temperature scale is not utilized on each figure. 



TABLE 5 



The number 
specified 


of periods at McMurdo Sound meeting each of the 
selection criteria. 


Selection Criteria 


Number of Days 
\Meeting Selection 
, Criteria 


Windspeed 

(mph) 


Cloud Cover 
(tenths) 


Length 
of Period 
(hours) 


L 5 


4. 5 


24 


68 


IN 

M 

O 


^ 5 


24 


156 


6. 15 


6. 5 


24 


244 


i 15 


f 5 


18 


357 


6 15 


£ 5 


36 


91 


15 


6. 5 


48 


50 


£ 15 


4 5 


60 


30 


± 15 


£ 5 


72 


22 


> 15 


> 5 


24 


78 


all 


4 5 


24 


710 


all 


> 5 


24 


458 


4 15 


all .. 


24 


607 


> 15 


all 


24 


183 
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TABLE 6 



The number of periods at the South Pole meeting each of 
the specified selection criteria. 



Selection Criteria 


Number of Days 
Meeting Selection 
Criteria 


Windspeed 

(mph) 


Cloud Cover 
(tenths) 


Length 
of Period 
(hours ) 


£ 10 






* 5 


24 


20 


£ 15 






£ 5 


24 


442 


CO 

o 






£ 5 


24 


828 


o 

CO 






£ 5 


24 


1217 


£ 45 






k 5 


24 


1217 


>15 






>5 


24 


86 


>30 






>5 


24 


0 


>30 






all 


24 


0 


>15, 


£ 


30 


£ 5 


24 


111 


>15, 


L 


30 


> 5 


24 


30 


> 20, 


c 


30 


> 5 


24 


2 


>20, 


£ 


30 


£ 5 


24 


2 
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7. Results of McMurdo Sound Data. 



Figure 6 shows the mean daily temperature variation de- 
termined by each observation time over all consecutive 0000 
to 2400 periods during the polar night. The range between 
the maximum and minimum values is 0.22 degrees (F) which 
leads to the authors' contention that there is no diurnal 
temperature variation at McMurdo Sound. The same curve com- 
puted by Barrigar with just six years of data had a varia- 
tion of 0.5 degrees (F) and it is felt that the mean diurnal 
variations reported by early investigators would not have 
been found, had additional data been available to them. Fig- 
ure 7 shows the mean temperature variation averaged at each 
observation time over all consecutive 48-hour periods. With 
only half as many intervals being averaged as in Figure 6, 
the temperature range has increased to 0.4 degrees (F). How- 
ever, there is no correspondence between the first and last 
24-hour periods of the curve as would be the case if a true 
diurnal temperature variation did exist. Note that the first 
2100 observation temperature is a relative maximum, while the 
temperature at the second is the curve’s minimum. 

Figure 8 represents the mean seasonal trend of the tem- 
perature, averaged daily using the ten years of data avail- 
able. The curve shows the general long term, downward trend 
expected during the polar night. The considerable aperiodic 
fluctuations reflect the variability of temperature during 
the polar night due to radiative as well as atmospheric cir- 
culation effects. 
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I Figure 6 

i 

Mean temperature variation at McMurdo Sound 
averaged for all days during the polar night 
for the years 1956-1965. Two temperature 
scales are illustrated. 




T 

E 

M 

P 

E 

R 

A 

T 

U 

R 

E 

(F) 



Figure 7 

Mean temperature variation at McMurdo Sound 
averaged for all consecutive 48-hour periods. 
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Figure 8 

Mean seasonal trend of the temperature at McMur 
Sound, averaged daily during the period 1956-19 
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Barrigar’s results supported Hisdal's suggestion that 
windspeeds ^ 15 mph and cloudiness L 5/10 were the optimum 
criteria for the selection of inversion type days. The re- 
sults of the authors' use of these selection criteria are 
shown in Figure 9. The curves have almost the exact same 
shapes as Barrigar's, although the new data increased the 
number of 24-hour periods meeting the criteria from 174 to 
244. The mean temperature curve exhibits a diurnal fluctua- 
tion and the Figure 9b statistical bias curve has the expected 
convex shape with an almost one degree range between its maxi- 
mum and minimum values. The Figure 9a meteorological tempera- 
ture curve does, however, show a diurnal variation of 1.1 de- 
grees (F) indicating the existence of a slight diurnal tem- 
perature variation for the days selected. 

Figures 10 and 11 show the results for selected days 
with windspeeds — 10 and - 5 mph respectively, and it is noted 
that as the selection criteria became more restrictive, the 
number of periods meeting the criteria decreased. It was 
anticipated that as the criteria for maximum allowable wind- 
speed was decreased, the inversion effect would become strong- 
er, resulting in increased concave curvature of the statisti- 
cal bias curves. The average temperatures of the mean tem- 
perature variation curves decreases as the windspeeds are re- 
duced, indicating intensification of the inversion, but the 
corresponding increase in the statistical bias effect is lack- 
ing as seen in studying Figures 9b, 10b and lib. This can 
most probably be attributed to the increased effect of random- 
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Figure 9a: Comparison of the mean 24-hour 

temperature variation with the column averaged 
meteorological and the statistical bias vari- 
tions at McMurdo Sound. 

(wind speed ^ 15 mph; cloudiness — 5/l0) 
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Figure 9b: Comparison of the mean 24-hour 

temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed — 15 mph; cloudiness £ 5/10) 
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Figure 10a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed £ 10 mph; cloudiness ^ 5/l0) 






Figure 10b: Conpc^rison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed ^10 mph; cloudiness ^5/ 10) 
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Figure Ha: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed —5 mph; cloudiness := 5/10) 
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Figure lib: Comparison of the mean 24-hour 
'temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMur'do Sound. 

(wind speed ^5 mph; cloudiness — 5/10) 
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ness in the data as the number of periods meeting the selec- 
tion criteria was reduced. 

The randomness effect is also evidenced in noting the 
increased variations between the statistical bias curves of 
the a and b figures. It is interesting to note the cases 
meeting the various selection criteria are few in number com- 
pared to the over 3600 possible 24-hour periods in the ten 
years of data. This number is based on the four starting 
times used in the study. 

In all three figures of this set, the mean temperature 
variation curves are quite similar and suggest a diurnal 
variation with a minimum temperature at 0600 and a maximum 
near midnight. The meteorological variation curves, however, 
do not show similar correspondence. Indeed they exhibit a 
variety of times of maximum and minimum temperatures, and 
cannot lead to the conclusion of a true diurnal temperature 
variation for inversion type days. 

The second set of McMurdo Sound curves, Figures 12 
through 16 result from the selection of periods of varying 
length having cloud cover - 5/10 and windspeeds - 15 mph. 
Periods of 18, 36, 48, 60 and 72 hours were selected to de- 
termine the applicability of the statistical bias effect over 
intervals other than the normal 24-hour day. 

From Figure 12b, it is seen that a smooth concave sta- 
tistical bias curve resulted for the 357, 18-hour periods 
that met the selection criteria. However, as the length of 
the interval is increased, and the number of periods meeting 
the criteria is decreased, the statistical bias curves become 
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Figure 12a: Comparison of the mean 18-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed ^ 15 mph; cloudiness - 5/l0) 



T 

E 

M 

P 

E 

R 

A 

T 

U 

R 

E 

(F) 




MEAN 

MET. 

BIAS 

357 cases 



HOUR OF OBSERVATION (LST) 



m*mniiHiai 



55 





Figure 12b: Comparison of the mean 18-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 
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Figure 13a: Comparison of the mean 36-hour 

temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed - 15 mph ; cloudiness - 5/l0) 
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Figure 13b: Comparison of the mean 36-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed —15 mph; cloudiness £ 5/l0) 
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Figure 14a: Comparison of the mean 48-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed - 15 mph ; cloudiness - 5/l0) 
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Figure 14b: Comparison of the mean 48-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed ^ 15 mph; cloudiness — 5/l0) 
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Figure 15a: Comparison of the mean 60-hour 
temperature variation with the column 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed - 15 mph ; cloudiness — 5/l0) 
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Figure 15b: Comparison of the mean 60-hour 
temperature variation with the diagonally 
averaged meteorological and the statisti- 
cal bias variations at McMurdo Sound. 

(wind speed — 15 mph; cloudiness -5/l0) 
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Figure 16a: Comparison of the mean 72-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed - 15 mph ; cloudiness ^ 5/l0) 
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Figure 16b: Comparison of the mean 72-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed £15 mph ; cloudiness ^ 5/l0) 
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more complex. For the 60 and 72 hour intervals, the bias 
curves have definite convex curvatures. It was expected 
that the bias curves would become less smooth as the number 
of qualifying periods decreased and the effects of random- 
ness in the data increased, but the convex curvature was 
unexpected and is not fully understood. 

An explanation may be that the periods of strongest 
inversion are of relatively short duration. This is not to 
say that the inversion is destroyed, but merely that its 
effectiveness is reduced by a random or cyclic increase in 
windspeed and/or cloudiness to values approaching the limit- 
ing values of the criteria. Increases above the criteria 
could, of course, occur between the times of observation. 

The effect of these increases may be sufficient to slow or 
even to temporarily reverse the otherwise longer term cool- 
ing process. A few such periods, randomly contained in the 
relatively few number of periods meeting the selection cri- 
teria for the longer intervals, could result in the observed 
bias curves. 

There is little variation in the average temperature of 
the observed temperature curves in this set of figures. The 
meteorological temperature variation curves for the 24 , 48 
and 72 hour intervals show a tendency for a relative minimum 
temperature to occur at 1500. There is otherwise, very little 
consistency in the times of occurrence of maximum and minimum 
temperatures in the meteorological variation curves of the 
set. Therefore, no support is given to the theory of a di- 
urnal temperature variation. 
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The next set of curves was computed to determine the 
applicability of a statistical bias effect for other than 
strong inversion periods and to provide the means for 
evaluation of the relative magnitudes of the effect of cloud 
and wind on the temperature inversion. It was anticipated 
that periods meeting the criteria of windspeeds > 15 mph and 
cloudiness > 5/l0 would be non- inversion in nature and would 
exhibit a convex statistical bias curve. Figure 17b, however, 
shows it to be concave during the first half of the period 
and rather flat for the remainder. A possible reason for this 
could be a lag between the onset of winds > 15 mph and cloud 
cover > 5/l0 and the destruction of the inversion. However, 
comparison of the mean temperature variation curve tempera- 
tures of Figures 9a and 17a show that in the present instance, 
the temperatures are 17 degrees warmer than, and probably 
thus physically removed from the inversion periods contained 
in Figure 9a. 

Figures 18 and 19 show the results for selection cri- 
teria of all cloud amounts, but with windspeeds L 15 mph and 
then > 15 mph respectively. With fewer cases the statistical 
bias curve of Figure 18a is not as smooth as that of Figure 
19b. The fact that the bias curve for windspeeds > 15 mph 
does not have convex curvature lends support to the sugges- 
tion that cloudiness, rather than windspeed plays the domi- 
nant role in the polar night temperature fluctuations. 

Figures 20 and 21 show the results for the selection 

criteria of all windspeeds, but with cloudiness k 5/ 10 and 
> 5/10 respectively. Comparison of the statistical bias 
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curves shows the marked change from concave to convex curva- 
ture as cloudiness criteria shifts from less than to greater 
than 5/l0 coverage. This curvature change is also accompanied 
by a 14 degree (F) warmer mean temperature curve and the con- 
siderable radiative warming effect of increased cloud cover 
is again illustrated. 

Four of the five meteorological temperature curves of 
this group of figures have a minimum temperatures at 0600 
and three have maximum temperatures at 1500. In none of the 
figures is the range between maximum and minimum temperatures 
more than 0.9 degrees (F) , and this is little support for 
implying existence of a diurnal temperature variation 
strongly suggested by the mean temperature variation curve of 
Figure 17a. 
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Figure 17a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed > 15 mph; cloudiness > 5/10) 
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Figure 17b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed > 15 mph; cloudiness > 5/10) 
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Figure 18a: Comparison of the mean 24-hour 
temperature variation with the column 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed > 15 mph; all cloud amounts) 
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Figure 18b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and statistical bias 
variations at McMurdo Sound. 

(wind speed > 15 mph ; all cloud amounts) 
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Figure 19a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(wind speed —15 mph; all cloud amounts) 
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Figure 19b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(wind speed -15 mph; all cloud amounts) 
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Figure 20a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(all wind speeds; cloudiness - 5/10) 
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Figure 20b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at McMurdo Sound. 

(all wind speeds; cloudiness £ 5/l0) 
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Figure 21a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(all wind speeds; cloudiness > 5/10) 
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Figure 21b: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at McMurdo Sound. 

(all wind speeds; cloudiness - 5/l0) 
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8. Results of South Pole Data. 



Figure 22 is the polar night mean daily temperature 
variation at the South Pole, averaged over all consecutive 
0000 to 2400 days for seven years. Two temperature scales 
are shown and the variation is just slightly more than two 
tenths of one degree (F). Only those days having all eight 
observations were averaged and because of the missing data, 

64 out of a possible 736 days were eliminated. The remain- 
ing 672 periods on which the curves are based are sufficient 
number on which to base the statement that there is no nor- 
mal diurnal temperature variation during the polar night at 
the South Pole. In light of the McMurdo Sound discussion, 
an increased number of cases should serve only to reduce the 
varience of the mean temperature variation curve. 

The mean daily seasonal temperature trend at the South 
Pole, Figure 23, showed even wider fluctuations than evi- 
denced at McMurdo Sound. The same considerable aperiodic 
fluctuations were noted on the computer print-outs of the 
seasonal trends for the individual years during this coldest 
time of year at the South Pole. The temperature fluctuations 
are thus quite active and it would be interesting to modify 
the computer programs to provide for a direct comparison of 
the mean windspeed, cloud and temperature variations. 

Examination of the data print-outs from the two stations 
showed that during the polar night the South Pole generally 
had much stronger winds than did McMurdo Sound. This is 
further evidenced in the first set of curves, Figures 24 
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Figure 23 

Mean temperature variation at the South Pole 
averaged for all days during the polar night 
for the years 1957-1964. Two temperature 
scales are illustrated. 
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Mean seasonal trend of the temperature at the 
South Pole, averaged daily during the period 1957-1964 
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through 27 for which were selected 24-hour periods with 
cloudiness — 5/10 and with windspeeds increasing from 10 
to 15 to 20 to 30 mph, respectively. Table 6 shows how the 
number of periods meeting the selection criteria increased 
from 20 to 1217 as the minimum windspeed requirement was 
raised. The few periods having windspeeds ^ 10 mph makes 
the representativeness of the curves of Figure 24 suspect, 
although smoothing indicates that the statistical bias curve 
is concave. The statistical bias curves of the remaining 
figures are markedly concave as expected and it appears that 
windspeeds 4 20 mph and cloudiness - 5/l0 provide optimum 
criteria for the selection of 24-hour inversion periods at 
the South Pole during the polar night. It is noted, however, 
that as the windspeeds are increased, the average values of 
the mean temperature curves also increase so that with cloudi- 
ness — 5/10, the periods of strongest inversion are those with 
the least amounts of wind. 

The mean temperature curves of Figures 25, 26, and 27 
are noted to be dominantly influenced by statistical bias 
curves. The meteorological temperature variation curves of 
these three figures all have relative minimums at 0300 and 
tend toward maximum temperatures at the beginning, middle, 
and ends of the period. However, since the periods selected 
accumulate from the succeeding figures and since the range 
between maximum and minimum temperatures is one-half degree 
(F) at most, there is no basis for assuming a diurnal tem- 
perature variation for these selected days. 
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Figure 24a: Comparison of the mean 24-hour 
temperature variation witji the column averaged 
meteorological and the statistical bias vari- 
ations at the South Pole. * 

(wind speed ^10 mph; cloudiness ^ 5/10) 
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Figure 24b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical bias 
variations at the South Pole. 

(wind speed — 10 mphj cloudiness — 5/10) 
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Figure 25a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at the South Pole. 

(wind speed ^15 mph; cloudiness -5/l0) 
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Figure 25b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at the South Pole. 

(wind speed ^ 15 mph; cloudiness ^ 5/l0) 
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Figure 26a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at the South Pole. 

(wind speed — 20; cloudiness £ 5/l0) 
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Figure 26b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at the South Pole. 

(wind speed — 20 mph; cloudiness — 5/l0) 
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Figure 27b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at the South Pole. 

(wind speed -30 mph ; cloudiness :=• 5/l0) 
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Figure 28 shows the results of selecting the 24-hour 
intervals with windspeeds > 15 mph and cloudiness > 5/10. 

The statistical bias curve is convex with a range of 2.5 
degrees (F). The number of periods selected has dropped to 
86, but the bias curve is quite smooth and typical of that 
expected for non- inversion periods. As noted previously the 
South Pole windspeeds were found to be higher than those at 
McMurdo Sound. But in trying various selection criteria it 
was found that with no restrictions on the cloud amount, 
there were no 24-hour periods with windspeeds consistantly 
greater than 30 mph. Additionally, there were no new 24- 
hour periods selected with cloudiness - 5/l0 when the wind- 
speed limit was increased from - 30 mph to — 45 mph. 

The last set of two curves, Figures 29 and 30 reveals 
the dominance of cloud over wind as a principal factor in 
warming of the surface layers at the South Pole. 24-hour 
periods with winds between 16 and 30 mph were selected for 
both figures, but in Figure 29, cloudiness was - 5/l0 and 
in Figure 29 > 5/l0. In Figure 29b, the statistical bias 
curve is nearly flat, but in Figure 30b it is concave and 
typical of the non- inversion situation. The mean tempera- 
tures were also slightly warmer in the latter non- inversion 
case. Windspeed remains a factor however, as can be noted 
from Figure 29 where the high windspeeds prevent the occur- 
rence of a general inversion situation. 
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Figure 28a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at the South Pole. 

(wind speed , >15 mph; cloudiness >5/10) 
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Figure 28b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at the South Pole. 

(wind speed > 15; cloudiness >5/10) 
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Figure 29a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at the South Pole. 

(wind speed > 15 and £ 30 mph; cloudiness — 5/l0) 
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Figure 29b: Comparison of the mean 24-hour 

temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at the South Pole. 

(wind speed > 15 and —30 mph; cloudiness — 5/l0) 




00 03 06 09 12 15 18 21 00 



HOUR OF OBSERVATION (LST) 



94 



Figure 30a: Comparison of the mean 24-hour 
temperature variation with the column averaged 
meteorological and the statistical bias vari- 
ations at the South Pole. 

(wind speed > 15 and —30 mph ; cloudiness "> 5 / 10) 
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Figure 30b: Comparison of the mean 24-hour 
temperature variation with the diagonally 
averaged meteorological and the statistical 
bias variations at the South Pole. 

(wind speed > 15 and — 30 mph; cloud iness> 5/10) 
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10. Conclusions and Acknowledgements. 



The averaging all of the available data did not reveal 
a diurnal temperature variation during the polar night at 
either McMurdo Sound or the South Pole. It is concluded 
that previously reported diurnal variations were due to 
randomness in the data as a result of the investigators not 
having available a sufficient amount of data. 

It was found that the mean temperature variations com- 
puted for specially selected days were considerably influ- 
enced by statistical bias, which when removed from the mean, 
left no consistent pattern of diurnal temperature variation. 
In addition, the data randomness is sufficient to warrant 
the use of both the column and diagonal averaging techniques 
in future studies. 

During the polar night, the marked fluctuations in the 
long term daily temperature trend seem to result more from 
cloud cover variation than from the more normal atmospheric 
circulation effects and it is recommended that further in- 
vestigation be conducted to determine more precisely the 
effects of cloud, windspeed and perhaps wind direction on 
the fluctuations. 

It is recommended that the computer programs be modi- 
fied to include the column and diagonal averaging computa- 
tions and also that the computer graphing sub-routines be 
employed to plot graphs. 

The authors wish to express their appreciation for the 
advice and guidance given during the conduct of this investi- 
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gation by Professor W. van der Bijl of the U. S. Naval 
Postgraduate School. 
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